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Abstract 

Albedo  is  the  ratio  of  total  hemispherical  reflected  (upwelling)  to  incoming 
(downwelling)  radiative  flux  (or  irradiance)  at  a  surface.  For  operations  on 
the  McMurdo  Ice  Shelf,  Antarctica,  albedo  is  a  controlling  factor.  Impuri¬ 
ties  such  as  dust,  soot,  mineral,  and  other  organic  deposits  on  a  snow  or 
ice  surface  can  dramatically  lower  albedo,  increase  solar  energy  absorption 
in  the  ice  shelf,  and  significantly  alter  the  energy  balance,  resulting  in  in¬ 
creased  melting,  snow  density  variations,  and  compromised  structural  in¬ 
tegrity  of  the  snow  and  ice  matrix.  The  occurrence  of  such  impurities  at 
Pegasus  Runway  may  have  been  a  factor  in  its  decline  and  replacement 
with  the  Phoenix  Runway  in  the  2016-  2017  field  season.  Therefore,  the 
National  Science  Foundation  (NSF)  requested  that  the  U.S.  Army  Cold  Re¬ 
gions  Research  and  Engineering  Laboratory  (CRREL)  establish  a  long¬ 
term  record  of  surface  albedo  at  the  Pegasus  and  Phoenix  Runway  sites.  To 
accomplish  this,  data  from  the  National  Aeronautics  and  Space  Admin¬ 
istration  (NASA)  Moderate  Resolution  Imaging  Spectroradiometer 
(MODIS)  satellite  was  used.  MODIS  satellite  data  has  a  daily  temporal  res¬ 
olution  and  a  significant  period  of  record  (2000-  present).  Additionally,  its 
narrowband  surface  reflectance  may  be  used  as  a  proxy  for  albedo.  This  re¬ 
port  documents  the  results  of  these  analyses. 
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1  Introduction 

1.1  Background 

McMurdo  is  the  main  transportation  hub  for  the  U.S.  Antarctic  Program 
(USAP)  on  the  Antarctic  Continent.  Passengers,  equipment,  and  supplies 
are  brought  into  McMurdo  by  vessel  and  air  and  then  redistributed  by  air 
and  overland  traverse  to  the  South  Pole  and  Antarctic  camps  across  the 
continent.  In  addition  to  a  dock,  McMurdo  maintains  airfields  to  support 
the  movement  of  personnel  and  supplies  into  and  out  of  Antarctica.  At 
McMurdo,  there  are  no  easily  accessible  locations  to  place  a  runway  di¬ 
rectly  on  soil;  as  a  consequence,  the  air  operations  into  and  out  of 
McMurdo  have  relied  on  establishing  and  maintaining  runways  (for 
wheeled  aircraft)  and  skiways  (for  ski  equipped  aircraft)  on  snow  and  ice. 

In  general,  the  integrity  of  the  Pegasus  and  Phoenix  Runways,  both  on  the 
McMurdo  Ice  Shelf,  depends  on  their  strength,  which  is  strongly  influ¬ 
enced  by  the  net  heat  transfer  into  the  runway  surface  during  the  brief  but 
intense  peak  of  austral  summer.  The  flux  of  downwelling  shortwave  solar 
energy  is  significant  during  the  austral  summer  when  24- hour- per- day 
sunshine  persists.  The  amount  of  this  radiative  flux  absorbed  by  the  run¬ 
ways  is  controlled  by  the  albedo  of  the  surface.  During  the  20 12-  20 13  and 
2013-  2014  seasons,  Pegasus  Runway  experienced  significant  reductions 
in  the  strength  of  the  snow  and  ice  layers  that  composed  the  runway,  lead¬ 
ing  to  extended  shutdowns.  In  addition,  during  short  warming  trends  that 
caused  temporary  reduction  in  the  strength  of  the  airfield,  operations  at 
the  Pegasus  Runway  were  shifted  to  the  "night"  when  the  sun  was  low  in 
the  sky  and  the  ice  was  slightly  colder  and  stronger.  Longer  warming 
trends  led  to  the  presence  of  liquid  water  on  the  runway  surface,  which 
significantly  lowered  the  runway  albedo.  The  formation  of  liquid  water  was 
a  tipping  point  for  runway  conditions.  After  this  point  was  reached,  the 
runway  degradation  was  rapid.  These  combined  factors  led  to  the  decline 
in  consistent  operations  at  Pegasus  Runway  and  ultimately  the  replace¬ 
ment  of  Pegasus  Runway  with  the  Phoenix  Runway  in  the  20 16-  20 17  field 
season  (Daly  et  al.  2015;  Haehnel  et  al.  2013). 

Impurities  such  as  dust,  soot,  mineral,  and  organic  deposits  on  a  snow  or 
ice  surface  also  lowers  the  surface  albedo,  increasing  solar  energy  absorp¬ 
tion.  This  increased  energy  uptake  can  initiate  melting  of  the  snow;  and  as 
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mentioned  above,  the  presence  of  melt  water  in  the  snow  further  reduces 
albedo  and  increases  heating  of  the  snow,  further  reducing  pavement 
strength  and  promoting  more  melting.  To  prevent  or  delay  reaching  the 
melt  tipping  point,  it  is  critical  to  maintain  the  albedo  of  the  snow  surface 
as  high  as  possible.  Based  on  some  qualitative  analysis  from  ground  obser¬ 
vations  and  contemporary  satellite  imagery  of  the  impacts  of  impurities,  a 
new  location  (referred  to  as  Phoenix)  was  identified  for  airfield  construc¬ 
tion  as  it  is  generally  further  away  from  a  zone  of  increased  dust  deposi¬ 
tion  and  closer  to  a  zone  of  increased  snow  accumulation  (Haehnel  et  al. 
2014). 

1.2  Objective 

Because  of  the  operational  difficulties  of  the  Pegasus  Runway  and  con¬ 
struction  of  the  new  Phoenix  Runway,  the  National  Science  Foundation 
(NSF)  requested  that  the  U.S.  Army  Cold  Regions  Research  and  Engineer¬ 
ing  Laboratory  (CRREL)  establish  a  long-term  record  of  surface  albedo  at 
the  two  runway  sites  (Figure  1) . 


Figure  1.  DigitalGlobe  WorldView-3  panchromatic  (grayscale)  satellite  collection  (25 
November  2015)  of  the  newly  established  Phoenix  and  existing  Pegasus  Runways.  Green 
and  red  polygons  represent  the  processing  window  for  analysis  (see  Section  2). 
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1.3  Approach 

Spectral  albedo,  a  (A),  is  defined  as  the  ratio  of  total  hemispherical  re¬ 
flected  (upwelling)  to  incoming  (downwelling)  radiative  flux  (or  irradi- 
ance)  at  a  surface  at  a  specific  wavelength.  A: 


a  (A)  = 


Fui  <U) 

Fa  (0,A) 


where 


Fu  =  the  upwelling  radiative  flux; 

Fd  =  the  downwelling  radiative  flux;  and 
0  =  the  surface. 

For  operations  on  the  McMurdo  Ice  Shelf,  Antarctica,  the  total  albedo,  or 
broadband  albedo,  at,  is  of  more  interest  as  it  is  a  measure  of  the  total  en¬ 
ergy  absorbed  by  the  ice  or  snow  surface.  This  is  found  in  terms  of  the 
spectral  albedo  as  (Perovich  1996) 

_  J  a(A)Fd(0,  A)dA 
Ut~  1  Fa  (0,  A)dA 

In  application,  this  equation  is  integrated  over  the  solar  (or  shortwave)  ra¬ 
diation  band  (approximately  300-  3000  nm)  to  estimate  the  broadband 
albedo. 

This  study  explores  using  the  National  Aeronautics  and  Space  Administra¬ 
tion's  (NASA)  Moderate  Resolution  Imaging  Spectroradiometer  (MODIS) 
satellite  data  to  determine  temporal  changes  in  surface  albedo  at  the  Pega¬ 
sus  and  Phoenix  Runway  sites.  MODIS  was  chosen  because  of  its  temporal 
resolution  (daily)  and  extensive  period  of  record  (2000-  present).  The 
MODIS  sensor  suite  exists  onboard  both  the  NASA  Terra  and  Aqua  satel¬ 
lites.  MODIS  Terra  data  are  presented  here  as  MODIS  Aqua  data  showed 
similar  results.  We  did  learn  that,  unfortunately,  the  MODIS  land- surface 
albedo  products  (MCD43A3)  are  not  available  because  the  region  of  inter¬ 
est  encompassing  the  Pegasus  and  Phoenix  Runways  has  been  masked  out 
(flagged  as  ocean)  from  the  processing  given  its  location  on  the  Ross  Ice 
Shelf,  which  floats  on  the  McMurdo  Sound. 
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Because  of  this  limitation,  MODIS's  narrowband  surface  reflectance 
(which  could  be  evaluated  as  individual  products  or  as  an  aggregate, 
broadband  reflectance)  was  used  as  a  proxy  for  albedo.  As  previously 
stated,  albedo  is  the  ratio  of  the  irradiance  reflected  from  the  surface  (leav¬ 
ing  at  all  angles)  to  the  total  irradiance  striking  the  surface  from  the  entire 
hemisphere  above  (incoming  radiation  at  all  angles)  at  all  wavelengths.  By 
contrast,  reflectance  is  the  ratio  of  reflected  radiation  observed  from  a  par¬ 
ticular  viewing  angle  (the  sensor  zenith)  and  illuminated  from  a  particular 
direction  (the  solar  azimuth  and  zenith) .  Additionally,  in  the  case  of  imag¬ 
ing  sensors,  reflectance  measurements  are  limited  to  bands  or  wavelength 
ranges.  Despite  these  differences,  reflectance  and  albedo  are  closely  re¬ 
lated.  Satellite- based  albedo  products  are  constructed  using  reflectance 
data  (Strahler  and  Muller  1999).  Because  collection  swaths  overlap  slightly 
and  shift  predictably  with  time  (MODI  S  Terra  cycles  through  its  orbital 
pattern  every  16  days),  each  day  there  are  several  observations  of  each 
ground  location,  each  with  a  different  sensor  and  solar  geometry  (and  dif¬ 
ferent  from  the  preceding  and  succeeding  day) .  By  aggregating  16  days  of 
data  and  interpolating  over  all  possible  angles,  end  member  albedo  cases 
("white  sky"  and  "black  sky")  can  be  modeled.  The  actual  albedo  or  "blue 
sky"  albedo  falls  somewhere  in  between  the  end  members  and  is  depend¬ 
ent  on  the  wavelength  range,  optical  depth,  aerosol  type,  and  solar  angle 
(Strahler  and  Muller  1999).  Because  albedo  is  derived  from  reflectance,  we 
assume  that  it  may  be  used  as  a  proxy  if  considered  across  several  wave¬ 
length  ranges  (bands),  particularly  over  longer  time  periods  ( >16  days) .  In 
general,  reflectance  values  should  follow  the  same  trend  as  albedo,  though 
the  magnitude  of  albedo  and  reflectance  may  differ. 
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2  Data  Set  and  Methodology 

The  raw  MODI  S  data  are  initially  collected  as  a  swath  product  (considered 
a  Level  1  product)  in  a  grid  that  is  aligned  and  unique  to  the  path  of  the 
sensor.  Higher- level  products  are  reprojected  to  a  "standard  grid,"  which  is 
very  distorted  at  high  latitudes.  For  example,  this  study  used  the  MODIS 
Surface  Reflectance  Product  (Level  2,  product  name:  MOD09GA),  which  is 
postprocessed  into  a  sinusoidal  projection  with  a  pixel  size  (at  the  equator) 
of  approximately  463  m  (NASA  LP  DAAC  20 16) .  At  the  extreme  latitude  of 
the  study  site  (77°51'  S),  there  is,  necessarily,  a  significant  difference  and 
distortion  in  the  layout  of  the  swath  data  and  that  of  the  Level  2  sinusoidal 
projection,  whose  ground  footprint  tapers  to  a  singularity  at  the  poles 
(NASA  LP  DAAC  20 16) .  Figure  2  illustrates  the  distortion  in  these  pixels 
by  the  comparison  to  WGS-84. 

Figure  2.  MODIS  pixels  (-500  m,  originally  in  sinusoidal  projection)  reprojected  and 
displayed  in  WGS-84  (geographic  projection)  are  illustrated  by  the  blue  grid { non¬ 
square  pixels)  overlaid  on  the  DigitalGlobe  Worldview-3  image  from  Fig.  1.  To 
minimize  distortion  and  errors  in  processing,  map  projection  uses  the  native  satellite 
collection  geometry.  Note  that  the  image  is  flipped  180°  when  compared  to  Fig.  1. 


To  mitigate  potential  error  caused  by  this  resampling  and  to  provide  a 
more  holistic  view  of  the  reflectance  over  each  of  the  Pegasus  and  Phoenix 
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Runways,  we  averaged  reflectance  over  the  pixels  that  intersected  the  dig¬ 
itized  extent  of  each  runway;  however,  because  of  the  orientation  of  the 
pixels  involved  (Figure  2),  a  simple  average  of  overlapping  pixels  is  not 
representative  of  each  runway's  reflectance  but  instead  of  a  much  larger 
area.  To  address  this,  we  calculated  the  percentage  of  each  pixel  that  actu¬ 
ally  overlapped  the  surveyed  extent  of  the  runway  (Figure  3)  and  used  that 
percentage  to  weigh  that  pixel's  contribution  to  the  runway  average: 

p  _  2n=i  wn^n 

—  v-’P 

In=i  Wn 

where 


Rr  =  the  percentage  reflectance  of  pixel  overlapping  runway; 

Rn  =  the  reflectance  of  pixel  n; 

wn  =  the  fraction  of  pixel  n  that  overlaps  the  runway;  and 
p  =  the  subset  of  pixels  for  each  individual  runway  that  are  doud- 
free,  high-quality  data.  Data  quality  is  determined  from  the 
MODLAND_  QA  "00"  flag  in  the  NASA  MODI  S  Terra 
(MOD09GA)  and  Aqua  (MYD09GA)  reflectance  products. 

Note  that  this  is  equivalent  to  a  simple  average  if,  for  all  pixels,  wn  =  1. 

Figure  3.  Raster  of  pixel  weights  (i vn),  or  the  percentage  of  an  overlapping  pixel  that 
is  inside  the  surveyed  extent  of  a  runway  (Phoenix  Runway  shown  center,  outlined  in 
black,  and  re-projected  in  the  M0DIS  sinusoidal  projection;  in  this  projection,  the 
M0DIS  pixels  appear  square). 
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3  Analysis 

The  surface  temperature  of  the  sun  generally  characterizes  the  solar  spec¬ 
trum.  This  spectrum  defines  the  spectral  irradiance  across  all  wavelengths 
of  radiation  as  the  Earth's  atmosphere  reduces  the  solar  irradiance  that 
reaches  the  Earth's  surface.  Ozone  (03),  water  vapor  (H2O),  and  carbon  di¬ 
oxide  (CO2)  absorb  this  radiation  at  certain  wavelengths.  Multispectral 
satellites  generally  leverage  regions  of  the  electromagnetic  spectrum  out¬ 
side  of  these  absorption  bands  (atmospheric  windows)  where  transmission 
is  greatest  through  the  atmosphere.  Figure  4  depicts  the  solar  irradiance 
outside  of  the  atmosphere,  the  direct  solar  irradiance  at  the  Earth's  sur¬ 
face,  and  the  NASA  MODIS  spectral  bands  that  typically  leverage  the  at¬ 
mospheric  windows  of  transmission.  This  study  analyzed  MODIS  bands  3, 
4, 1,  and  2  (listed  in  order  of  increasing  wavelength)  as  narrowband  and 
broadband  analogues  to  albedo  based  on  the  highest  relative  amounts  of 
solar  irradiance. 


Figure  4.  Solar  spectral  irradiance  across  the  electromagnetic  spectrum  outside  of  the 
atmosphere  (green)  and  at  the  Earth’s  surface  (red).  Discrete  NASA  MODIS  satellite 
bands  are  highlighted  in  orange. 


For  available  MODIS  Terra  satellite  observations.  Figures  5-  8  highlight 
the  average  daily  reflectance  for  the  period  of  record  (2000-  2016)  by  us- 
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ing  the  processing  methodology  for  calculating  representative  pixel  per¬ 
centages  as  explained  in  the  previous  section.  The  figures  show  for  the  fol¬ 
lowing  spectral  ranges  the  results  of  pixels  encompassing  both  the  Phoenix 
and  Pegasus  runways:  459-  479  nm  (MODIS  Band  3),  545-  565  nm 
(MODIS  Band  4),  620-  670  nm  (MODIS  Band  1),  and  841-  876  nm 
(MODIS  Band  2). 


In  general,  for  each  MODIS  band,  the  pixels  encompassing  the  Phoenix 
Runway  exhibit  a  higher  surface  reflectance  than  the  pixels  encompassing 
the  Pegasus  Runway.  This  is  analogous  to  higher  surface  albedo  across  the 
broadband  wavelengths.  Moreover,  these  data  are  consistent  with  the  solar 
spectral  irradiance  curves  in  Figure  4  such  that  shorter  wavelengths  show 
relatively  higher  spectral  reflectance  than  longer  wavelengths  for  both  run¬ 
ways  based  on  the  corresponding  higher  solar  spectral  irradiance  at 
shorter  wavelengths. 


Figure  5.  Average  daily  reflectance  over  2000-2016  for  pixels  encompassing  Phoenix  ( red) 
and  Pegasus  ( blue )  Runways  for  MODIS  Band  3  (459-479  nm). 
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Figure  6.  Average  daily  reflectance  over  2000-2016  for  pixels  encompassing  Phoenix  (red) 
and  Pegasus  (blue)  Runways  for  MODIS  Band  4  (545-565  nm). 
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Figure  7.  Average  daily  reflectance  over  2000-2016  for  pixels  encompassing  Phoenix  (red) 
and  Pegasus  (blue)  Runways  for  MODIS  Band  1  (620-670  nm). 
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Figure  8.  Average  daily  reflectance  over  2000-2016  for  pixels  encompassing  Phoenix  (red) 
and  Pegasus  (blue)  Runways  for  M0DIS  Band  2  (841-876  nm). 
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4  Conclusions 


We  found  that  over  the  period  of  observation  (2000-  2016),  the  MODIS 
data  showed  that  generally  the  broadband  reflectance,  and  by  analogy  al¬ 
bedo,  is  higher  at  the  Phoenix  site  than  at  the  Pegasus  Runway.  While  it  is 
difficult  to  directly  compare  these  satellite  observations  to  the  range  of  al¬ 
bedo  from  fresh  Antarctic  snow  and  bare  blue  or  sea  ice,  these  results  are 
reasonable  when  compared  to  other  general  studies  of  snow  and  ice  albedo 
as  well  as  Antarctic  specific  studies  (Gardner  and  Sharp  2010;  Wuttke  et 
al.  2006).  The  MODIS  data  shows  obvious  seasonal  variations  in  the  al¬ 
bedo  with  the  highest  albedo  observed  in  the  spring  (September),  a  steady 
decline  through  the  summer,  and  a  rebound  in  late  March  and  early  April 
(Figure  9) .  This  is  to  be  expected  as  summer  snowmelt  and  metamorpho¬ 
sis  of  the  snow  pack  will  reduce  the  albedo  of  snow. 


Figure  9.  Monthly  average  reflectance  over  2000-2016  for  pixels  encompassing  Phoenix 
(red)  and  Pegasus  (blue)  Runways  for  MODIS  Band  3  (459-479  nm). 
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In  Figure  10,  the  MODIS  data  also  show  that  both  sites  have  shown  a 
steady,  but  slow,  downward  annual  trend  in  albedo  from  2000  to  2012. 
There  appears  to  be  a  more  rapid  decline  in  albedo  at  Pegasus  from  20 13 
to  20 16.  A  similar,  but  less  pronounced  decline  in  albedo  may  also  be  oc¬ 
curring  at  the  Phoenix  site  likewise  starting  in  about  20 13.  Additional 
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analysis  is  required  to  assess  the  significance  of  these  trends  and  how  they 
may  persist  in  future  operational  seasons. 

Figure  10.  Annual  average  reflectance  over  2000-2016  for  pixels  encompassing  Phoenix 
(red)  and  Pegasus  (blue)  Runways  for  MODIS  Band  3  (459-479  nm). 


Annual  Average  (2000-2016)  Reflectance  (MODIS  Terra) 

Band  3  (459-479  nm) 

- PEGASUS  - PHOEIMIX - Linear  (PEGASUS) - Linear  (PHOENIX) 


Additional  work  may  also  include  ( 1)  further  investigation  of  how  well 
broadband  MODIS  reflectance  bands  estimate  trends  in  snow  albedo  and 
(2)  whether  changes  in  albedo  may  be  leveraged  toward  estimating  when 
the  melt  tipping  point  of  the  runway  is  reached  at  each  site.  Challenges  ex¬ 
ist  overcoming  the  coarse  spatial  resolution  of  the  satellite  (MODIS  ~  500 
meters)  and  optical  methods  being  impacted  by  doud- cover.  This  results 
in  a  noisy  signal  (both  spatially  and  temporally)  and  requires  a  statistical 
approach  performed  over  a  time  series  (seasonal  to  decadal) .  Additional, 
consistent  in- situ  observations  of  albedo  and  melt  conditions  are  required 
to  determine  whether  this  type  of  satellite  remote  sensing  technique  can 
link  lowered  surface  reflectance  and  surface  melt  toward  a  runway  mainte¬ 
nance  quality  control  tool. 
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